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ABSTRACT 



Aims. We examine the interactions of various instabilities in rotating stars, which usually are considered as independent. 
Methods. An analytical study of the problem is performed, account is given to radiative losses, /i-gradients and horizontal 
turbulence. 

Results. The diffusion coefficient for an ensemble of instabilities is not given by the sum of the specific coefficients for 
each instability, but by the solution of a general equation. We find that thermohaline mixing is possible in low-mass red 
giants only if the horizontal turbulence is very weak. In rotating stars the Rayleigh-Taylor and the shear instabilities 
need simultaneous treating. This has for consequence that rotation laws of the form l/r a are predicted to be unstable 
for a > 1.6568, while the usual Rayleigh criterion predicts instability only for a > 2. Also, the shear instabilities are 
somehow reduced in Main Sequence stars by the effect of the Rayleigh-Taylor criterion. Various instability criteria 
should be expressed differently in rotating stars than in simplified geometries. 

Key words, physical processes: mixing, stars: structure, stars: evolution, stars: rotation 



1. Introduction 

Various instabilities and processes intervene in stars, partic- 
ularly in rotating stars, and produce some mixing in radia- 
tive zones. These processes are often expres sed by diffusion 
coefficients (lEndal fc Sofia!lll7l[Zahl[l992h . Currently the 
global diffusion coefficient is expressed as the sum of the 
particular coefficients of the various effects considered in- 
dependently. However, there are many physical interactions 
between the different instabilities, with possible amplifica- 
tion or damping effects. Our aim here is to consider the 
interaction of different instabilities and their associated cri- 
teria. 

The effects and criteria we consider here are the ther - 
mohaline instability (|Ulrichlll972T:lKippenhahn et al.lll980l) . 
the shear instability expressed b y the Richardson crite rion, 
the Rayleigh-Taylor instability (|Chandrasekharl[l96ll ). the 
semiconvective instability, the instability charact erized by 
the Solberg-H0iland criterion ( Wasiutv nskil 1 194(1 ) and the 
baroclinic insta bility expressed by the Golre i ch-Schubert - 
Fricke criterion (jGoldreich fc Schubertl Il967l: iFrickel 1 1968(1 . 
Account needs to be given also to non-adiabatic effects, 
to the //-gradients and to the horizontal turbulence gener- 
ally present in differentially rotating stars. The interactions 
of horizontal turbulence with shear instabilit ies and merid- 
ional c irculati on were studied respective ly bv I Talon fc Zab.nl 
(|1997ft and bv lChabover fc Zahnl (|1992ft . The exact ampli- 
tude of the horizontal turbulenc e is uncertain and various 
estimates haye bee n performed (|Zahnl[l992t iMaederl feOOS: 
iMathis et al"1l2004l) . 
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An interesting consequence of this work is that some 
of the well known instabilities should never be considered 
alone in rotating stars, but always in their interactions with 
other ones. This leads to modified expressions of some well 
known criteria when applied to rotating stars. 

In Sect. 2, we start by deriving the basic expression 
of thermohaline mixing in two ways, both showing that a 
proper account of the mean molecular weight is often miss- 
ing in many derivations. We also examine analytically the 
effect of the horizontal turbulence on thermohaline mixing. 
In Sect. 3, we show that many instabilities are interacting 
and should not be treated separately, but are better anal- 
ysed in the frame of a more general equation. In Sect. 4, we 
reconsider some well known stability criteria and provide 
new expressions for them in rotating stars. Sect. 5 gives 
the conclusions. 

2. The double - diffusive mixing revisited 

2.1. The thermohaline mixing 

We first consider the effects of the gradients of mean molec- 
ular weight ii on the stability, since such gradients are al- 
ways present in stellar interiors as a result of nuclear reac- 
tions and ionization. Most of the time the /Lt-gradients have 
a stabilizing effect, since generally regions with higher fi are 
deeper in the star. 

The /x-gradient may also be destabilizing, when in a 
medium subject to gravity some upper layers have a higher 
mean molecular weight /i than deeper layers (V/i < 0). 
The medium may nevertheless be stable with respect to 
the Ledoux criterion for convection. This is possible, if the 
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upper layers with higher fi are hot enough to allow a sta- 
ble density gradient. However, if heat diffusion is fast (faster 
than matter diffusion) , the medium of higher /i has the time 
to cool, become denser and start sinking. This is the ther- 
mohaline instability, which takes the form of " salt fingers" 
in lab experiments. 

The interest for thermohaline mixing has been revived 
over recent years in relation with the observed chemi- 
cal abundances of red giants, which require some signifi- 
cant internal mixing at the level of the H-burning shell. 
ICharbonnel fc Zahni (|2007l ) and ICharbonnel fc Laeardel 
(|2010f) investigated the effects of thermohaline instability 
in red giant stars using a classical treatment for the asso- 
ciated diffusion coefficient; they found in particular an ex- 
cellent agreement between theoretical predictions and the 
observations of Li, C, and N isotope s on the upper red giant 
branch. Also, lLagarde et al.l (|2012t ) showed that the inclu- 
sion of thermohaline mixing in stellar models provides a 
very elegant solution to the long-stan ding helium-3 prob- 
lem in the Galaxy. On the other side, ICantiello fc Langerl 
(|2010H and IWachlin etall (|2011h showed that to explain 
the surface abundances observed at the bump in the lu- 
minosity function, the speed of the thermohaline mixing 
process needs to be more some orders of magnitude higher 
than in their models. These different results are due, at 
least for a large part, to the different geometries adopted 
for the fluid elements. In these last two works, the fluid 
elements are essentially spher ical, while oth er authors like 
ICharbonnel fc Lagardd (120101) in the line of lUlrichl (l!972h 
consider elongated "salt fingers" . We will see below (Eq. [3] 
and [5]) that the ratio of the vertical length of the salt fingers 
to their diameter intervene to the square as a multiplicative 
factor in the thermohaline diffusion coefficient. 

The theory has been analyzed by 
iDenissenkov fc Pinsonneault! (|2008F) in relation with 
3 He burning in the envelopes of red giants, a work con - 
firmed by 2D numerical simula tions (IDenissenkovl l2010f). 
Furth er numerical simulations (Deniss enkov fc Merrvfieldl 
1 2 llh suggest that the thermohaline mixing is not efficient 
enough (by a factor of 50) to account for the observed com- 
position anomalies of r ed giants. Numerica l 3D simulations 
have been performed (jTraxler et al.l 1201 if) , these authors 
develop new prescriptions and support the inability of 
thermohaline mixin g to achieve the observed mixing, while 
I Brown et al.l (|2012T) find more transport than suggested 
by Traxler et al. The effect of radiative levitation and 
tur bulence on thermohaline mixing are nicely analyzed 
by IVauclair fc Theadd (|2012h . On the whole, we see a 
large variety of results, often related to the choice of some 
physical parameters. 



2.2. The account of the chemically inhomogeneous 
surroundings 

Often a perturbation method is used to der ive the coeffi- 
cient of thermohaline diffus ion in the line of Uhich (|l972j ) 
or a blob method following iKippenhahn et al.l (jl980ll . The 
problem is that the possible ambient variations of /i are 
generally neglecte d . Let us consider the blob method. 
Kip penhahn et afl (|l980l ) have shown that the falling ve- 
locity Vfj, of a blob with an excess A/i > with respect to 
the surrounding medium is given by 
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where r t hcrm is the thermal adjustment timescale given by 
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where K is the thermal diffusivity K = 4 a cT 3 /(3 k Q 2 Cp) 
and d the diameter of the blob. The factor of 6 applies to 
a spherical bubble, while elongated blobs like salt fingers 
have larger values, e.g. 12, 20 or more. The usual notations 
are used, in particular Hp the pressure scale height, 6 = 
— (din g/dlnT)p^ , <p = (91ng/91n/i) PT . In a radiative 
medium V a d > V, thus for Afj, > one has < and the 
blob is sinking. 

The thermohaline diffusion coefficient is Z?thi = 
(because the transport is anisotropic, there is no factor 
1/3), thus one has 
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There, for the expression of the relative excess of mean 
molecular weight, the following expression is generally 
A/i//i = I (dlnfi/dr) = —IVn/Hp, where I is the mean 
free path of the fluid element. The geometrical factors 
6 (i/d) 2 can be grouped in a constant Cthi, which may be- 
come large for salt fingers. A n expression of the form ([3]) has 
been found bv lUlricbT (|1972l ) with a coefficient Cthi = 685, 
a similar expression but with Cthi = 12 was obtained by 
IKippenhahn et al.l ljl980T ). The sa me expression but with 
Cthi = 1000 has been used by ICharbonnel fc Lagardd 
(|2010f ). implying that the mixing-length of the fluid ele- 
ments is very large with respect their diameters. Similar 
expressions, but wit h different coefficients, wer e also ob- 



tained and used by IDenissenkov fc Pinsonneault (2008), 
iTheado fc Vauclairl (|2012h and lVauclair fc Theadol (|2012h . 

A critical point is that expression (JTJ) and consequently 
([3]) apply for thermohaline mixing in a surrounding homoge- 
neous medium. This is not the general case, e.g. in a radia- 
tive shell H-burning region, there is a /i-gradient through- 
out the zone where mixing develops. (We may also remark 
that in the derivation of ((3j) it is a bit strange to assume the 
surrounding medium homogenous and at the same time to 
make use of a certain V/x. Maybe, this V/i, could be con- 
sidered as the average gradient between the blob and the 
surrounding medium, but this is not very satisfactory). 

Thus, let us consider an inhomogenous medium with a 
certain /x- gradient. the g eneral form of the sinking velocity 
thus (|Maederll2009r) . 



v u is 
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As a consequence, the thermohaline diffusion coefficient is 

(~Vm) 
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This expres s ion supports the consistent finding by 
IDenissenkovl (|2010h of a term V/x at the denominator of 
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(O, a term which is generally ignored, both in theoretical 
papers and modeling of stellar evolution. We examine the 
size of the various terms in the numerical model of a 1.25 
Mq of solar composition, with an initial rotation velocity 
of 110 km s _1 at the stage of the re d giant bump with a lu- 
minos ity of 105 L®, as in Fig. 9 by ICharbonnel fe Lagardd 
(2010). The zone where thermohaline mixing is present ex- 
tends over about 2.6 % of the radius around a mass fraction 
Mr/M = 0.30 at the base of the extended outer convective 
zone. Over most of this small zone, (V a d— V) is around 0.15, 
while the term (<p/6) is about 3 • 10~ 6 , meaning that, 
even if V M is the leading term at the numerator its effect 
at the denominator is negligible here. As a result the coef- 
ficient of thermohaline diffusion is high, lying between 10 7 
and 10 10 cm 2 s _1 , while the coefficients resulting from other 
diffusion sources are in the range of 10 2 to 10 3 cm 2 s _1 . So 
far so good, this would confirm the importance of the ther- 
mohaline mixing, however the thermohaline instability, like 
other instabilities, is also influenced by non-adiabatic ef- 
fects and turbulence, which we now consider. 

2.3. The effects of thermal losses and turbulence 



the coefficient of horizontal turbulence D n plays the same 
role as the thermal diffusivity in the previous expression (|7f 
and write, 



, int 
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With proper account of heat losses and horizontal turbu- 
lence, the modified Richardson criterion for shear instability 
becomes (jTalon fc Zahnlll997D . 
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with A/ a 2 d = |L(V ad -V), and TV 2 = gv„ . (11) 



The s tandard value of lZi c — 1/4. Following Tal on fc Zahnl 
(|!997l ). we write with x = v£/6, 
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Let us first consider the effects of radiative thermal losses. 
In a medium with gravity, the net acceleration (by unit of 
length) due to the density difference between a fluid element 
and the surroundings is, 



d(S In g) 
dr 
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where the two parentheses on the right contain the dif- 
ferences between the internal and external gradients. The 
difference of T-gradients depends on th e radiative los ses 
determined by the thermal diffusivity K (|Maederlll995l ). 



Vint - V = 



T (V ad -V) and r=^, (7) 



where T is the ratio of the energy transported to the 
energy lost on the way (i.e. T equals the Peclet number 
over 6), with v and £ the typical velocity and mean free 
path of the motions. 

Let us now turn to turbulence, which originates mainly 
from differential rotation, both vertical (producing vertical 
shears) and horizontal (producing horizontal shears creat- 
ing horizontal turbulence). 

The (vertical) shear instability is produced by a sig- 
nificant difference Av of velocity between nearby layers of 
different vertical coordinate r. In a shear, there is an excess 
energy by unit of mass and surface, 



The heat transport by D n adds its effects to K . Assuming 
that the shear diffusion coefficient £> s hcar = (1/3)?; £ = 2x 
(the transport concerns t he vertical dir e ction) is small 
with respect to K and Z? n , iTalon fc Zahnl ([l997) have ob- 
tained Dshear with account of horizontal turbulence and 
heat losses. Let us note that this applies to a spherical fluid 
element. If not, the result for x has to be multiplied by 
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2.4. Limit for thermohaline instability 

If there is no energy available from the differential rotation 
for shear instability, the second member of (TT2"|) is zero. 
Then, if we also ignore the horizontal turbulence, we just 
get the limit 
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This is the same as ((3]), it also shows the presence of a 
term depending on the /j,-gradient at the denominator. 
For thermohaline mixing, V/x and N 2 are negative, thus 
the diffusion coefficient is evidently positive. As mentioned 
above, N 2 at the denominator increases the mixing. 

If some horizontal turbulence is present, but no signifi- 
cant shears, we get from (|T2|) 



(8) Ahi = 6 



locally present in the differential motions with respect to 
the case of an average velocity. The Richardson criterion 
defines the occurrence of shear instability over the zone in r 
considered. If the excess energy in the shear is greater than 
the work necessary for overcoming the stable temperature 
and /i-gradients, we get the usual Richardson criterion (see 
also Sect. S3]). 

As to the horizo ntal turbulence , it int ervenes by chang- 
ing the ^-gradients. ITalon fc Za hn (1997]) have shown that 



-KNl 



N 2 



(14) 



The first term in the parenthesis is positive and the second 
negative. Thus, for mixing to occur, one must have a limit 
on the value of the horizontal turbulence Z?h, 



D h < K 



-N 2 
I* 



N 2 



(15) 



A similar expression, but without the term in 
A^ 2 at the denominator, has been obtained by 
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iDenissenkov fc Pinsonneaultl (|2008l ). In general, the 
numerical calculations show that D n is smaller than K. 
However, for thermohaline mixing to occur, it is necessary 
that at the beginning (before the slow descending motions) 
the system is stable, thus N 2 d » \N?\, this means that 
the r.h.s. of (|T5]> may be much smaller than K. Thus, -D n 
becomes a factor which may cancel the large parenthesis 
in (|14l) and stop the mixing. This is consistent: when the 
differences of /i are small with respect to the stabilizing 
effect of the T- gradient, turbulence may easily kill the 
thermohaline mixing. 

As a numerical example, we consider the 1.25 Mq model 
discussed in Sect. 12.21 In this example, the r.h.s of (IT51) is 
1-2 orders of magnitude larger than Z? n , implying that the 
adopted horizontal turbulence does not inhibit the thermo- 
haline mixing. However, the situation may not be so sim- 
ple. As mentioned, there is a great uncertainty in the size 
of the horizontal tur bulence and three different expressions 
have b een proposed (IZahnlll992|lMae der 2003; lMathis et aTl 
120041) . ICharbonnel fc Lagardd (|2010D have taken the lowest 
bound of the first one. The last two, likely more physical 
(|Maederll2009h . give similar numerical values, larger by 4 
to 5 orders of magnitude than the first one. If applicable, 
such high values of -D n are sufficient to kill the thermoha- 
line mixing. Thus, the problem of the thermohaline mixing 
remains open and we see that a key to its solution lies in 
a better knowledge of the horizontal turbulence in stars, in 
addition to the proper choice of the geometrical parameters 
of the fluid elements. 



2.5. Semiconvective diffusion 

We now turn to semiconvection, for which the expression 
of the diffusion coefficient is related to that of the thermo- 
haline mixing. Formally, semiconvection occurs when the 
thermal gradient V is intermediate between the condition 
of stability predicted by the Ledoux criterion and the in- 
stability predicted by Schwarzschild's criterion, 



< V < V i; 



(16) 



Vi n t being very close to V ac j in stellar interiors. The ve- 
locity of the semiconvective mixing is also given by (|4]) 
as for thermohaline mixing. The blobs being in dynami- 
cal equilibrium, the equation of state implies ^f- = j ^ , 
where AT for a blob displaced over a distance Ar is 
Ml = ( Vint - V) With account of © and Q a coef- 
ficien t of semiconvective diffusion may be obtained ([Maederl 

MM), 



2T 
f+T 



K{V- V ad ) 
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For r tending towards infi nity (adiabat i c case ), one gets 
the coefficient obtained by lLanger et al.l (|1983l ). However, 
semiconvection arises from nonadiabatic effects and in stel- 
lar interiors, T lies between 10~ 2 and 10~ 3 . Thus, the nu- 
merical coefficient in front of (TT7| is of this order of a 
magnitude. 




Fig. 1. Geometry data concerning the Rayleigh- Taylor in- 
stability. 



3. Instabilities in rotating stars 

3.1. The coupling of the Rayleigh-Taylor and shear 
instabilities 

The effects of shears, thermal diffusion, thermohaline mix- 
ing and horizontal turbulence are all already included in 
Eq. (TT2"j) . However, we must also account for the stabiliz- 
ing or destabilizing effects related to the distribution of the 
specific angular momentum j = vo 2 f2. Indeed, if a (de- 
scending) salt finger brings more specific angular momen- 
tum inwards than present in the surrounding medium, this 
difference of j will tend to bring the fluid element back to 
its initial position after some oscillations of frequency Nq, 
thus the medium is stable. At the opposite, if the salt finger 
brings less specific angular momentum than present in the 
surroundings, it tends to continue its way and an instabil- 
ity is produced. An analogous discussion can be made for 
outwards motions. This is well known from the Rayleigh 
instability criterion: a medium with j decreasing outwards 
is unstable. 

In order to examine closely this effect, let us consider a 
star in differential rotation. For now, we need not to specify 
more the dependence of the angular velocity fl on the spa- 
tial variables. Let us consider a fluid element at a distance 
wq and let us move it outwards, as illustrated in Fig. [TJ 
to a very close point at distance w. The difference of the 
centrifugal forces acting on the blob by unit of mass may 
be written 
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The work <57^ e ntrif delivered by the bubble to the medium 
during this small displacement is 

J_ d(il 2 uj 4 ) 
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duj 



(zu — ZUqY 



(19) 



If fl decreases very fast with increasing w, may be neg- 
ative. This means that the motion of the bubble sponta- 
neously delivers some work 57^ C ntrif to the medium, which 
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is thus unstable. At the opposite, if Nq is positive, one has 
to give some energy to the bubble in order to move it out- 
wards. The medium is stable. This leads to the Rayleigh 
criterion for the stability of a distribution of angular mo- 
mentum, 



i d(n 2 w 4 ) 



dm 



> 



(20) 



We emphasize that the above effects of the work made 
by the centrifugal force leading to Rayleigh criterion (|2"0"|) 
and the effects of shears as expressed by ([8]) are two 
different physical effects. In some particular hydrodynamic 
situations, these two effects can be treated separately. 
For example, in plane parallel flows with different ve- 
locities, the problem of centrifugal force is not relevant. 
At the opposite, in a differentially rotating star, both 
shear effects and angular momentum effects may influ- 
ence the stability and should then be treated simultaneously. 

The displacement of a fluid element implies the ex- 
change of bubbles, leading to some mixed fluid element with 
an average spatial velocity. As recalled above, this makes 
an energy (1/4) (Sv) 2 available by unit of mass. In a star 
with differential rotation of the "shellular" form O(r), the 
velocity difference 5v is 



ov = r smtf — — or 
dr 



(21) 



where d is the colatitude of the point considered (cf. Fig. 
[I]). Thus, an energy ^7^hcar is made available during an 
exchange, 



<5T sh ca r = Tr 2 sin 2 tf(^) 5r 2 
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If, this energy available from the shear is smaller than the 
energy (I20p necessary to destabilize the angular momentum 
distribution, the medium is stable, i.e. if 
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In the equatorial regions for a case with Sl(r), this can be 
written, 



2 — — ] < 4 N? 



'dsi\ 

or in terms of SI 



(24) 



If SI decreases outwards faster than this limit, the medium 
is unstable. Rotation laws of the form l/r a with a be- 
tween 1.6568 and 2.0 are predicted to be stable by the usual 
Rayleigh criterion, while if we also account for the excess of 
energy available in the shear motions they are unstable. It 
is clear that there is no reason to ignore the excess of energy 
available in shear motions, which contribute to destabilize 
the system. Therefore in a differentially rotating star, the 
Rayleigh criterion must be modified (see Sect. I4.1j) . 

Reciprocally, when shear instability is considered in a 
rotating star, the stabilizing or destabilizing effect of the re- 
distribution must also be accounted for. Below, the two ef- 
fects are coupled and considered simultaneously. This leads 
to some modifications of well-known criteria. 

3.2. A more general coupling 

Let us try to be more general and consider simultane- 
ously various effects which may influence the stability of 
the medium in a rotating star with Sl(r): 

— the effects of thermal gradients, 

— the thermohaline mixing, 

— the semiconvective diffusion, 

— the shear mixing to the local excess of energy in differ- 
ential rotating layers, 

— the stabilizing or destabilizing effect of the distribution 
of angular momentum, 

— the radiative losses, 

— the transport of heat by the horizontal turbulence, 

— the element diffusion in the medium due to the horizon- 
tal turbulence. 

Following Eq. dTUl) and the previous subsection, also taking 
into account that centrigugal force operates orthogonally 
to the rotation axis, we write for the complete form of the 
instability expression, 
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N. 



ad 
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where the standard value of lZi c is 1/4 and we have ac- 
counted for the effect of the distribution of the specific 
angular momentum j on the stability. Making the same 
transformation x = v£/6 as in (fl"2"]) . we get the following 
instability condition, 



dsiy ^ 

i- I -3- J < 16 ft 2 



irfi 



dSl 
dr 



(25) 



The equality of the two members of this expression gives a 
quadratic equation in (dSl/dr), the solution of which is 



dSl\ _ n 
dr J r 



(i±V5) 



(26) 



To be consistent with a solution in absence of shear, we 
take the sign minus in the above expression. This means 
that in order to be stable, fi(r) should not decrease faster 
than a limiting ^-distribution, 
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This general equation may be used to get the total resulting 
diffusion coefficient D to t = 2x of the chemical elements in 
a star, taking into account the many effects involved. The 
stability condition may be written as a quadratic equation 
of the form Ax 2 + Bx + C > 0, 



Nf, 



with 



N 2 (K + D h )+N 2 _ Sv (K + 2D h ) 



N 2 _ Sv (D h K + D 2 ) 



> 0, 



(30) 
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i d(n 2 w 4 ) 

U7 3 dw 



sin z9 — lZi c 



(31) 



Expression of N^_ Sv is a modified form of the Rayleigh 
oscillation frequency, accounting for both the angular mo- 
mentum distribution and the excess of energy in the shear. 
Eqn. (|3T)|) is the general equation, which should generally 
be considered in a rotating star to account for the vari- 
ous effects mentioned above. It contains the Schwarzschild, 
the Ledoux, the Solberg-H0iland, the Rayleigh-Taylor, 
the Richardson and the GSF criteria with radiative losses 
and horizontal turbulence all in one single expression. 
Simplifications are possible, but they should be justified. 

In the first evolutionary phases, is generally likely to 
have a stabilizing influence since the fi-gradient is not very 
steep, this means that the instabilities could be somehow 
reduced in the equatorial regions. In the advanced phases 
after the end of the He-burning phases, N% may be desta - 
bilizing over very limited regions dHirschi fc Maederll2010lK 
there equatorial instabilities may be favoured. 

At the equator, coordinates w and r coincide, the effect 
of the centrifugal force is maximum, while at the pole this 
effect disappears and vanishes. The distribution of the 
instability inside the star is no longer spherically symmet- 
ric. Indeed, all properties of rotating stars are two dimen- 
sional: the shape, the isobars, the meridional circulation 
and the various instabilities have a dependence in colati- 
tude t9. The usual practice (jKippenhahn fc Thomasl[T970f ) 
is to consider the equations of structure at an average ra- 
dius rp, such that the volume of the sphere of radius rp is 
equal to the real volume of the oblate rotating star. This 
is also v alid for differentially rot ating star with shellular 
rotation (M evnet fc Maedenll997[ ). The radius rp is equiv- 
alent to the radius of the oblate star at colatitude d given 
by A(cos?9) = 0, i.e. cos 2 = 1/3, which corresponds to 
i9 = 54.7 degrees. Thus, in Eq. d3TJ) , we put sin?? = y/{2/3) 
and we may thus keep the advantages of the ID modeling. 
However, all these procedures may satisfactorily apply for 
the lower rotation velocities, while for high velocities, there 
may be significant differences with the complex 2D reality 
of rotating stars. 



3.3. The general solution 

In stellar models, it is often considered that the total dif- 
fusion coefficient Aot is the sum of several diffusion coeffi- 
cients A, each one corresponding to an instability i, 



Act = 2 A 



(32) 



Such a summation of coefficients A is undoubtedly wrong! 
The solution of an equation of the second degree like 
is not the sum of the limits of some particular cases. The 
solution is that of the complete equation. This has for con- 
sequence to introduce correlations between the various ef- 
fects in the form of products of terms appearing in the co- 
efficients A, B and C of the quadratic equation (j30|) . These 
are 



A=[K 



(33) 



B = N^ d D h +N%(K + A) + N*_ 5v (K + 2D h ) , 
C = N&_ Sv (D h K + Dl). 



The solution of ([3"U| brings coupling of all intervening fac- 
tors. We notice in particular the cross-products of the heat 
losses and of the horizontal turbulence with the effects of 
the thermal, fx— and fi-gradients. There are also products 
of these three gradients and also direct products of heat 
losses and horizontal turbulence. This indicates that the 
full solution of the diffusion coefficient D — \2x\ depends 
on the coupling of the various physical effects considered 
and is by no means a sum like 



4. Some stability criteria reconsidered 

The above results lead us to reconsider some basic stabil- 
ity criteria in differentially rotating star, which indeed only 
account for one specific aspect of the instability and ignore 
other ones, unavoidably present in rotating stars. In par- 
ticular, we see that the complete equation (I3U|) contains 
the term Nq_ Sv . In agreement with the result of (|3.ip . the 
term Nq_ 5v should always be considered instead of Nq in 
the relevant criteria. 



4.1. A modified Rayleigh-Taylor criterion 

The usual Rayleigh-Taylor criterion for stability of a rotat- 
ing flow (in absence of gravity) is 



iV 2 > 0,with iV 2 



— 



i d(n 2 w A ) 

~~3 



dvo 



(34) 



It expresses that in order a rotating medium is stable, the 
angular momentum must increase with the distance to the 
rotation axis. 

If in the general equation (1501) we put K and A equal 
to zero and if we also do not consider the effects of thermal 
and /x-gradients, we are left with the stability condition 

Nh-sv > ■ (35) 
At a colatitude $ in a star, we should consider this new form 
of the criterion for the stability of differential rotation, If 
positive, it means that the excess energy in the differen- 
tial rotation is unable to destabilize a stable distribution of 
angular momentum. If negative, the medium is unstable. 

Indeed, the Rayleigh-Taylor criterion (f34|) may predict 
stability of a given differential rotation law, while in fact it 
would be unstable according to ([33)1 as discussed in Sect. 
13.11 The usual Rayleigh-Taylor criterion (l34l) should never 
apply as such in a rotating star, since the shears also favour 
instability. 

4.2. A modified form of the Solberg-H0iland criterion 

Let us first recall the usual Solberg-H0iland criterion for 
convective stability, 

NL+tf + NZ sintf > 0. 



(36) 

This to some extent combines the Ledoux and the Rayleigh- 
Taylor criteria. It expresses the stability with respect to 
convective motions (at the dynamical timescale B? / (GM)) 
with also account for the stability of the rotation law. 

If in the general equation Q30p we put K and A equal 
to zero, the terms B and C in (|34D are zero and we are left 
with the condition 



A 



N 2 aA + Nl + Nf 



n-5v 



> 



(37) 
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or explicitely 



4.4. The GSF criterion 



N. 



ad 



9 1 d(n 

N -| v 



dm 



sin I? 



i dv 
Hi* I 3- 

dr 



> 0. 



(38) 



This the modified form of the Solberg-H0iland criterion for 
stability, accounting for the fact that both the distribution 
of j and the energy excess in differential rotation are con- 
sidered simultaneously. The last term on the left hand side 
may reduce stability. Thus in a differentially rotating star, 
a situation considered as stable according to (j3"6")l could in 
reality be unstable. 

The Solberg-H0iland criterion does not account for ra- 
diative losses, which would introduce the term T/(T + 1) in 
front of the first term in (|3"6")l and (|3"5|) . This would bring 
a dependence on the thermal timescale. It neither account 
for a viscosity due to turbulence, which would introduce a 
dependence on the viscous timescale. The Solberg-H0iland 
instability essentially obeys the dynamical timescale. The 
terms accounting for radiative losses and (turbulent) vis- 
cosity will appear below in the so-called GSF criterion. In 
such a case, one would have a kind of semiconvective diffu- 
sion like (|17|) . including in addition the effects of rotation. 
The characteristic timescales is determined by the thermal 
timescale rather than by the dynamical timescale. 

4.3. The Richardson criterion 

The usual expression of the Richardson criterion for shear 
instability in a star with shellular rotation is 



V;; R /< | ^ 

dr 



(39) 



It expresses that if the excess energy in the shear overcomes 
the restoring energy available in the stable temperature- 
and /i-gradients, then mixing occurs. 

If we ignore thermal losses and horizontal turbulence, 
the modified form of the Richardson criterion is identical to 
the above modified expression (|3"5|) of the Solberg-H0iland 
criterion, since the shear effects have to be considered to- 
gether with those regarding the stability of the angular mo- 
mentum distribution. 

If account is given to the thermal losses with the term K 
and not to the horizontal turbulence, the stability relation 
(|3"0|) becomes, 



K 2 d 

[NlK- 



Nl 



x ■ 



n-8v\ 
K] > 0. 



(40) 



The corresponding diffusion coefficient -D s hcar = 2ir is, 



The Goldreich-Schubert-F ri cke ( GSF) criterion 
(|Goldreich fc Schubert! Il967t iFrickd I1968D considers 
the instability of a stellar rotation law subject to thermal 
diffusivity K and to viscosity v. Instability occurs for each 
of the two conditions, 



K 



< 



and 



an 2 

T 

dz 



> — Nr 



K 



T,ad ' 



(42) 
(43) 



The first inequality expresses the same idea as the Solberg- 
H0iland criterion, with account of thermal diffusivity. The 
second relation expresses the instability related to the dif- 
ferential rotation in the direction parallel to the rotation 
axis. 

The ratio v/K can also be expressed in term of V as in 
Sect. 12.31 where we see a term T/(r + 1) in front of A^ d . 
The viscosity is v = (1/3) v£ and T = v£/(6K). Formally 
one has the correspondence T = v/(2K) , the difference 
between a term like in (|42|) and terms like T/(r + I) comes 
from two facts: 

-I. In expressing T, account is given to the geometrical 
shape of the fluid element supposed to experience radiative 
losses, while in (|4"2"j) this i s not made. 

- 2. Also, as shown bv lMaederl ( | 19951 12009D . in the gen- 
eral case of radiative losses one must consider the ratio 
r/(r + 1) in front of A^ d . It is only for an extreme conduc- 
tivity (r — > 0) of the medium that one may just write V or 
v/K. At the opposite, in the perfect adiabatic case, r — > oo 
and the term in front of A^ d is 1. 

From these two remarks, one sees that the following 
writing of the standard GSF criterion is more appropriate, 



and 



r + i. 

dz 



m M + K sinz? < 

r 



> 



r + i 



(44) 
(45) 



with the same notations as in Sect. 12.31 

The effect of /^-gradient can also be accounted for, as 
well as the diffusivity of chemical elements, due in partic- 
ular to the horizontal turbulence. The instability is then 
triple diffusive. In a rotating star, the horizontal turbulence 
is a major diffusive proce ss, its effects on the GS F insta- 
bility have been studied (|Hirschi fc Maederll2010[) . As for 
previous criteria, we also need here to account for the ex- 
cess energy present in the shears. We do it in the same way 
as before in replacing Nq by Nf l _ Sv . The first expression of 
the GSF instability criterion is thus 



D, 



shear 



= 2K 



N 2 



N 2 



(41) 



This is the diffusion coefficient for the effects of shear, with 
account of the T- and /i-gradients, of the thermal losses 
and for the Rayleigh-Taylor instability. The complete form 
of the stability criterion, including thermal losses and hor- 
izontal turbulence, is given by (1301) . 



r 



N, 



r + i 

The second is 

on 2 



T. ad 



> 



Nf, + N, 



^T,ad 



n-8i 



< 



Nf, 



(46) 



;(§)<«> 
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These conditions are easier to realize than the standard 
ones, which means that the domain of the GSF instabil- 
ity is somehow extended by the accoun t of the excess en- 
ergy in the shear. As shown previously (jHirschi fc Maeder] 
2010), the GSF instability can be locally very large in pre- 
supernova stages, however the zones over which the GSF 
instability is acting are very narrow. It may be interesting 
to examine whether the above expressions are changing the 
results. As a matter of fact, we see that the first criterion 
is equivalent to our general stability equation (|50)l . The so- 
lution is given by the root of this quadratic equation. 

5. Conclusions 

We have considered simultaneously various instabilities in 
differentially rotating stars and given a general equation 
expressing the instability conditions with respect to an 
ensemble of effects, including horizontal turbulence. This 
introduces some coupling between the instabilities. We 
suggest that in rotating stars some instabilities should 
never be dissociated. In particular the instability of 
rotation law expressed by Rayleigh criterion and the shear 
instability should always be considered simultaneously. 
Modeling the instabilities in 1 D models of stars with low 
or moderate rotation velocities is possible by considering 
the effects at the root of P2(cost9) = 0. In some future 
work now in preparation, we will examine the consequences 
of the present approach on the mixing in stellar interiors. 

Acknowledgements: G.M., N.L. and C.C. acknowledge dif- 
ferent supports from the Swiss National Fund. 
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